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X , Abstract. We study the effects of Supernova (SN) feedback on the for- 

^ ■ mation of disc galaxies. For that purpose we run simulations using the 

I extended version of the code GADGET-2 which includes a treatment of 

■ chemical and energy feedback by SN explosions. We found that our model 

succeeds in setting a self-regulated star formation process since an impor- 
tant fraction of the cold gas from the center of the haloes is efficiently 
heated up and transported outwards. The impact of SN feedback on 
' galactic systems is also found to depend on virial mass: smaller systems 

are more strongly affected with star formation histories in which sev- 
eral starbursts can develop. Our implementation of SN feedback is also 
successful in producing violent outflows of chemical enriched material. 

Resumen. En este trabajo se estudia el impacto de explosiones de Su- 
pernova (SN) en la formacion de galaxias de disco. Para ello se utilizan 
simulaciones numericas realizadas con una version del codigo GADGET-2 
que incluye los efectos de enriquecimiento quimico e inyeccion de en- 
ergia debido a explosiones de SN. Este modelo es capaz de generar un 
mecanismo auto-regulado de formacion estelar, como resultado del ca- 
lentamiento eficiente del gas en el centre de los sistemas y su trans- 
porte hacia el halo. En nuestro modelo, el impacto de las explosiones 
de SN en sistemas galacticos depende de la masa virial de los mismos: los 
mas pequerios son afectados mas fuertemente pudiendo presentar varios 
brotes de formacion estelar. Nuestro modelo es tambien capaz de pro- 
ducir fuertes flujos de gas hacia las regiones externas de las galaxias, los 
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cuales transportan una importante fraccion de metales hacia el medio 
intergalactico. 



1. Introduction 

Supernova explosions are thought to play a fundamental role in the evolution 
of galaxies. Through the ejection of metals and energy into the interstellar 
medium, they significantly affect the gas collapse, the star formation activity 
and the chemical patterns. As a consequence, SN feedback would be able to set 
a self-regulated mechanism for the star formation activity and the enrichment 
of the Universe. Moreover, SN driven outflows are believed to be responsible for 
the presence of metals detected in the intergalactic medium (e.g. Tytler et al. 
1995; Songaila & Cowie 1996). 

In the context of numerical simulations, the treatment of SN feedback has found 
severe complications. Several works have developed models of different complex- 
ity to treat SN feedback in Smoothed Particle Hydrodynamics (SPH; Gingold Sz 
Monaghan 1977; Lucy 1977) codes (e.g. Katz & Gunn 1991; Navarro & White 
1993; Metzler k Evrard 1994; Marri & White 2003; Springel k Hernquist 2003), 
finding different levels of success. 

2. Numerical Models 

In this study, we analise simulations of isolated disc galaxies of different virial 
masses from idealized initial conditions. The initial conditions correspond to 
spherical distributions of superposed gaseous and dark matter particles, per- 
turbed to give rise to a density profile of the form p(r) oc r~^. The spheres are 
initially in solid body rotation with a spin parameter of A 0.1. We have used 
an initial number of 9000 both gaseous and dark matter particles for all masses. 
The simulations were performed with the extended version of the Tree-PM SPH 
code GADGET-2 (Springel & Hernquist 2002) which includes a treatment of 
chemical enrichment and energy feedback by Supernova (Scannapieco et al. 
2005a, b,c). The model is tight to a multiphase scheme which allows a better 
description of diffuse gas in the context of the SPH technique. The interested 
reader is referred to Scannapieco et al. 2005a,b for details on the numerical 
scheme. All experiments presented in this study include chemical enrichment 
and metal-dependent cooling, and they all have been run in the context of the 
multiphase treatment for the gas component. 

It is relevant to this study to mention that the energy feedback model has only 
one free input parameter. Within the model, we define for each star particle 
with associated SN explosions two gaseous phases in the following manner: a 
cold phase defined as the gas with T < 2T* and p > O.lp* {T^, = 4 x 10^ K, 
p* = 7 X 10~^^ g cm~^) and a hot phase formed by the remaining gas. The 
energy released in the explosions, as well as the metal content, are distributed 
into the cold and hot phases: the cold phase receives a fraction while the hot 
phase gets the remaining e/j = l—Cc- Hence, the parameter €c (or equivalently) 
is the only free parameter to assume. Note that the value of Cc is related to the 
power of the feedback. However, we note that the effects of varying its value are 
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not completely linear, owing to the non trivial interplay among energy release, 
heating and cooling (see Scannapieco et al. 2005b). 

3. Results 

We analise here results for isolated disc galaxy simulations, in the virial mass 
range 10^'^ — 10^^ Mq. In Fig. 1 we show the evolution of the star formation 
rate for tests of 10^^ Mq (upper panel), lO^'^'^ Mq (middle panel) and 
Iq9.5 ^-1 jYj^ systems (lower panel). We have plotted the curves for experiments 
performed without including energy feedback (solid lines) and including feedback 
with an input parameter of = 0.5 (dashed lines). These trends indicate that 
the inclusion of the model of SN energy feedback succeeds in regulating the star 
formation activity in galaxies. This is the result of the effects of the injection 
of energy into the interstellar gas: when SNc explosions take place, the gas is 
efficiently heated up and accelerated outwards. As a result, the cold gas density 
from which stars are formed decreases and the star formation rate is reduced. 
Note also that lower mass systems are more strongly affected, as expected since 
their potential wells are less efficient in retaining baryons (Larson 1974). 
As it was mentioned before, SN feedback is responsible of heating up the gas 
in star forming regions and, consequently, of generating important outflows of 
material. Since these outflows originate in the center of the systems where the 
gas is highly enriched with metals, they can transport an important fraction of 
the chemical elements into the haloes and even beyond. This could in principle 
explain the observed presence of metals in the intergalactic medium. As an 
example, we show in Fig. 2 the edge-on projected gas metallicity distribution 
for our test of 10^^ Mq virial mass system just after the starburst. From this 
figure we can see that our model is able to enrich the regions outside the centre 
of galaxies. Also note that the generated outflows are mostly perpendicular to 
the disc plane. 



4. Conclusions 

We have presented results of the effects of SN feedback on the evolution of 
isolated disc galaxies of different mass, by using an extended version of the 
code GADGET-2 which includes a treatment of chemical and energy feedback 
by Supernova. We have shown that the model succeeds in generating a self- 
regulated star formation process, as a consequence of the triggering of strong 
outflows of material. 

We found that the impact of SN feedback on galaxy formation depends on virial 
mass. Smaller systems are more strongly affected since they have shallower 
potential wells and hence arc less efficient in retaining baryons. In particular, 
we found that the star formation rates of small systems such as dwarf galaxies 
can show a series of starbursts, in consistency with recent observational findings 
(Kauffmann et al. 2003; Tolstoy ct al. 2003). 

In our model, galactic outflows transport an important fraction of the metal 
content from the centre to the outskirts of galaxies. This is of crucial importance 
since it is observed that the intergalactic medium is contaminated with heavy 
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Figure 1. Evolution of the star formation rates for our tests of iso- 
lated rotating spheres with virial masses of 10^^ Mq (upper panel), 
jQio.5 ^-1 (middle panel) and 10^-^ h"^ (lower panel). We have 
plotted the curves for experiments performed without including energy 
feedback (solid lines) and including feedback with an input parameter 
of ec = 0.5 (dashed lines). 



Supernova Feedback in SPH simulations 



5 



T = 0.73 Gyrh"' 




-60 -40 -20 20 40 60 

X [ /j ' kpc J 



0.01 0.10 1.00 




Figure 2. Edge-on gaseous metallicity distribution for our 1012 h-^ 
Mq virial mass system run including the feedback model with an input 
parameter of ec = 0.5 (dashed lines), after 0.75 Gyr of evolution. 
The metallicity scale is also shown. 
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metals. SN driven outflows such as those generated by our model constitute a 
natural explanation for this observational result. 
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